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ABSTRACT 

The purpose of this paper is to suggest a technique 
of cluster analysis which is similar in aim to the Interactive 
Intercolunuiar Correlation Analysis (IICA) , though different in 
detail. Two methods are proposed for extracting a single bipolar 
factor (a "contrast coropenent'J) directly from the initial 
similarities matrix. The advantages of this general approach are 
that: (a) It helps avoid certain misclassification problems inherent 

in IICA; (b) It is related in a straightforward way to conventional 
techniques of multidimensional scaling and therefore allows a unified 
treatment of dimensional and "typal" structures ; and (c) It provides 
an interesting solution to the problem of relations among linear 
contrasts based on different subsets of the stimuli. (Author/DJ) 
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1. Page 9, equation (26) 
now reads: 

"F(b) = b'lfe - [tr(R)/n] (Vl)" 
should read: 

”F(b) = b'Rb - [tr(R)/n] (b'l)^" 



2. Page 9> paragraph immediately following equation (26'), 
second sentence 
now reads: 

"This solution is invariant vinder transformations 
of the form r^j -*■ ar^^j + c , while the Method A 
solution is invariant only under transformations 
of the form r^j -► •" 

shovild read: 



"Solutions by either Method A or 
invariant under transformations 
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Michael Shafto 
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Introduction 

L. L. McQuitty (I 967 ) has suggested a technique of hierarchical cluster 
analysis called Iterative Intercolumnen* Correlational Analysis (IICA). McQuitty 
and Clark have provided a discussion of the mathematics of this technique, its 
application to real and artificial data, and its advantages and disadvantages 
in comparison with other methods (Clark & McQuitty, 1970} ^fcQuitty, 197'1} 

McQuitty, Abeles, & Clark, 1970; McQuitty & Clark, 1968 ). Coles and Stone 
( 1972 ) have suggested a related technique. 

IICA begins with a raw data matrix from which a first-order similarities 
matrix R^^^ is computed (or, in some cases, R^^^ niay be obtained directly 
by subjects' similarity, judgments). The larger the ij -th entry in V, 
the more "alike" or "similat," in some sense, stimuli i and are judged 
to be. A second-order similarities matrix R^^:^ is then computed by intercor- 
relating the columns of R^^^ . Ihat is, the ij -th entry in R^^^ is the 
product -moment correlation between columns ^ and of R^ ^ . Then R^ ^ 

(2) 

is computed by intercorrelating the columns of R' y and so on, until a 
matrix is produced in which all elements have absolute value unity, 

Trathin reasonable tolerance limits. The stim\ili are then partitioned into 

( k) 

tvfo subsets: Those with +1 in the first column of R' go in one subset; 

those with -1 go in the other. ( Any column coiild be used, not just the first.) 
A discussion of the convergence problem may be found in Clark and McQuitty (1970). 

^Research supported by National Science Foundation Grant GB 8025X with 
Princeton University, project on "Mathematical Techniques in Psychology," Harold 
Gulliksen, Principal Investigator, and by Educational Testing Service. 

The author is grateful to Dr. Walter Kristof for many valuable suggestions 
and comments bn an earlier version of this report. The assistance of Mrs. Ann^ 
King in supervising the preparation of the manuscript, and of Mr. Terry Birch in 
drawing the figures, is also gratefully acknowledged. 
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The purpose of this paper is to suggest a techniciue of cluster analysis 
which is similar in aim to IICA, though different in detail. Two methods will 
be proposed for extracting a single bipolar factor (a "contrast component") 
directly from the initial similen:ities matrix. The advantages of this general 
approach are that (a) it helps avoid certain misclassification problems 
inherent in IICAj (b) it is related in a strai^tforwen:d way to conventional 
techniques of multidimensional scaling (Torgerson, 1958) therefore allows 
a unified treatment of dimensional and "typal" structures, and (c) it provides 
an interesting solution to the problem of relations among linear contrasts 
based on different subsets of the stimuli. This last problem was initiaiLly 
raised by McQuitty ( 1967 ). 

I. Matrix Algebra of One IICA Iteration 

The following discussion of the matrix algebra of one IICA cycle is 
intended to clarify the relationship between IICA and the new techniques out- 
].ined belov;. 

Consider a typical iteration, starting with ^ and ending with ^ • 

The superscript will be dropped from R^^^ for purposes of the follo\^ing dis- 
cussion. 

R is an n X n symmetric, matrix which is positive semidefinite for 
k > 2 . Therefore, there exist matrices U and D , such that D is diagonal, 
U’ = U"^ , and R = UDU’ . 

Let 1 be a colximn vector of n 1 *s, and let I be the n x n identity 
matrix. Define 
(1) M = I - ll’/n ; 



( 2 ) 



R = MR 



( 5 ) 

(M 

( 5 ) 



A = R’R = UDU’M'MUDU* = UDU'MUDU' , since M is idempotent; 
T = U'MJ ; and 
b = U’l . 



Thus, 

(6) T = I - ^'/P* ^ 

and the general element of T is 

0) ■ "iV" ’ 

where 8. . is the Kronecker delta. 

10 

By the Cauchy-Schwarz inequality, 

(8) jb^l <'/n , i = l,...,n 



Therefore, 

(9) =l,...,n . 

Nov 7 suppose that the inequality in (8 ) holds for each i , as it almost 
always will with real data. Then all the diagonal elements of T are positive 
Define 



( 10 ) 

( 11 ) 

Thus, 

( 12 ) 



C = Cdiag(T)]“2T[diag(T)]"2 ; and 

‘ JL 

W = UDCdiag(T)]2 



A = WCW’ , 



and the 



k + 1 -order similarities matrix is given by 
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(15) = [diag(A)]"^A[diag(A)] 



Nov 7 C is positive semidefinite, since (4) and (lO) imply 




where M is idempotent, therefore positive semidefinite, and, by the "law of 
inertia," the transformations from M to C preserve definiteness. Further- 
more . it is clear from (10 ) that each diagonal element of C is unity • Thus , 
the necessary and sufficient conditions for C to be a matrix of cosines between 
pairs of vectors in Euclidean n -space are satisfied. 

Therefore, (12) is the familiar expression for the inner-products matrix j 



of a set of vectors , where the columns of W represent the coordinates of the 



In effect, then, the IICA method performs a transformation of the vectors 
(or points) that represent the judged stimuli. The nature of this transforma- 
tion is as follovfs : The vectors at "time k " had coordinates UD with respect 

to n orthogonal axes. The vectors at "time k + 1 " have coordinates 
[diag(A)]"% with respect to n oblique axes. The cosines between pairs of 
oblique axes are the elements of C . 

The effect of this transforriiation can be seen more clearly by noting 



vectors with respect to n oblique axes • • • ^31] 

cosine of the angle between w. and w . 

1 J 



,w ) , and c . . is the 
-n'' ' io 




Jl - b?/n , and 



(16) 






i / d 
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In equation (15 ), since = u^l , the expression under the radical 

attains its maximum precisely when tne sum of the elements of is 0 . 

Now d^^ = > which is the squeu:ed length of the vector u^^ . (If R 

v/ere the dispersion matrix of a set of n random variables, then this quantity 
would be the variance of the particular linear combination of those variables 
represented by u^ .) Tlius, in order for the elements of w^ to be "large," 
the length of u^ must be "large," and the sum of the elements of u^^ , i.e., 
the sum of the projections of the vectors representing the judged stimuli ^ 
u. , must be "small." 

^ 2 

Consider equation (l6); c. . is indeterminate O/O if either b^ - n 

or = n but this will seldom occ\ar with real data. As b^ ot b^ , 

u 

approaches 0 , so does c^^ . If Ib^l and Ib^l are large, and if they 
have the same sign, then c. . becomes large and negative j if they have 

X j 

opposite signs, then c. . becomes large and positive. 

Intuitively, equations (15) and (l6) represent two "processes" which are 
being "applied" simultaneously to the vectors which represent the judged stimuli 
in IICA. Equation (15) states that bipolar axes are "lengthened" while nonbi- 
polar axes are "shrunk." Equation (l6) states that bipolar axes tend to remain 
orthogonal to one another, while nonbipolar axes are rotated toward or away 
from one another so that they tend to "collapse" into bipolar axes. 

This is hov 7 IlCA converges toward a single bipolar axis , as illustrated 
in Clark and McQuitty (l970). In the early iterations the stimulus -vectors are 
transformed toward bipolarity. The greater the initial departure from bipolar- 
ity, the greater the "correction factors." As bipolarity is attained in the 



^Note that, if b^ = n , then bj^ = 0 , for all k / j 
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later liberations, the IICA process becomes sijnilar to the well-known power- 
method (Anderson, I958) for extracting the largest latent root and correspond- 
ing latent vector of a real symmetric matrix. The "process" represented by 
(16) becomes negligible in the later iterations. 

II. Alternative Methods 

The fact that IICA transforms the stimulus -vectors themselves, rather than 
providing a solution in terms of the original configuration, seems, a priori , 
to be a drawback. Are the clusters revealed by the method prominent in the 
data, or are they "weak"- -perhaps even artifacts of the method itself.' Besides 
such theoretical questions, there are practical problems with IICA, as shown in 
Section IV below, which can be avoided by the alternative techniques suggested 
here. Moreover, these alternative techniques allow the unified treatment of 
dimensional and "typal" structures, as originally suggested by McQuitty (1967 )• 

Two methods will be proposed. Neither of these methods requires additional 
assumptions about the initial data matrix. Both involve extracting a single 
bipolar factor in such a way as to display bhe major clusters, and both can be 
applied recursively to yield hierarchical solutions. Neither makes any trans- 
formation of the original stimulus configuration. 

Method A 

Ignoring equation (16), and concentrating on equation (l 5 )> we seek a 
vector b , such that b’Rb is "large" and Ib'lj is "small." (Note that 
this b is not the b of Section I. ) Proceeding rather directly, we seek to 
maximize b'Rb under the constraints b'^ = 1 and ^*1 = 0 . Introducing 
Lagrange multipliers 7^ , 7 g ^ we write 



O 
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(17) F(b,7i,72) = - 7i(b'b - 1) - 72(b’l) v 

Differentiating F with respect to b , ^2 ’ setting the 

derivative equal to 0 , yields 

( 18 ) 2Rb - 27 ^b - 72i = 0 5 

(19) b'b -1=0 ; and 

(20) b'l = 0 . 

Eremultiplying (l8) by 1' gives 

(21) 72 = 21'R^n . 

Premultiplying (l8) by b* gives 

(22) 7^ = b*Rb 
Substituting for 7 g in (l8) gives 

(25) (I - ll7n)Rb = 7^b , 

or, following our previous definition of M , 

(24) MRb = 7^b . 

But b'l = 0 by (20), and it is easy to show that b'l = 0 iff Nto = b . 

p 

Therefore, (24) is equivalent to 



(25) MR^to = 7^b 

From (19), (22), and (25), it follows that the desired solution, b* , is 
the normalized latent vector of MRM corresponding to the largest latent root. 
But what is MRM ? It is simply the scalar-products matrix of the stimulus - 



^This step, which shortens the derivation of the solution by about one page, 
was suggested by Dr. Walter Kristof. 
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vectors with respect to a coordinate system with origin at their centroid (see 
Torgerson, 1958^ P* 258, equation l4). 

The vector b* represents the most sadient linear contrast between sub- 
sets of the stimuli* The stimuli can be ordered with respect to the correspond- 
ing elements of b* • Two criteria whicii could then be used to partition the 
stimuli into subsets (clusters) are 

1. VJeak Criterion ; Examine the n - 1 differences between adjacent 
elements of b* (having arranged these elements in order of magnitude )> find 
the largest such difference, and make the split between the corresponding 
stimuli. This should suffice when the clusters are fairly distinct. 

2. Strong Criterion ; Consider each of the n - 1 possible splits 
between pairs of stimuli which are adjacent with respect to b* • For each 
possible split, the original similarities matrix R can be partitioned into 
similarities within clusters and similarities between clusters. Thus, for 
each split, a quantity can be computed which reflects the adequacy of the 
partition. It is naturally desirable to have large similarities within 
clusters and small similarities between clusters. Therefore, one formula 
which could be used would be formally identical to the formula for the altJia 
level by a median test. Choose the split which minimizes the "alpha" for the 
appropriate "one-tailed test." Of course, it is not suggested that the minimum 
"alpha" reflects the statistical significance of the clustering. It simply pro- 
vides an intuitively appealing objective function which is sensitive to cluster 
size as well as to differences in magnitude of similarities. 

Certainly other methods of partitioning could be suggested. Examination 
of b* itself, however, will often indicate the presence or absence of clear 
clusters, and should provide a check on the adequacy of any method of partitioning 
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Basically, what is bein.'j suggested is to use a one -dimensional scaling solution 
as a heuristic to reduce the number of possible partitions to be considered. 

* 

Method B 

The following method constrains lb’ll to be "small"--but not necessarily 
0 . Since we want to make b'Rb "large" and Ib’lj "small" at the same time, a 
natural function to maximize is simply b’Rb - (b'l) , again under the con- 

straint b 'b = 1 . But since the largest possible value of b'Rb is the 
largest latent root of R , which can be no larger than tr(R) , and the 
largest possible value of (b'l)^ is n , a more "balanced" objective function 
is 

(26) F(b) = b'Rb - [tr(R)/n] (b'l) 

= b'(R - [tr(R)/n]ll’)b , 

or, letting R* = R - [tr(R)/n]^' , 

(26’) b'R*b , b'b =1 . 

The dtisired solution is simply the normalized latent vector corresponding 
to the largest latent root of R*- . This solution is invariant under transforma- 
tions of the form r. . ar., + c , while the Method A solution is invariant only 

under transformations of the form r^^ -► ar^^ . 

The extraction of two or more vectors of R* or MRM may be useful when 
cross-classification, rathei than general hierarchical clustering, is desired. 
Cross-classification is a special case of general hierarchical clustering, 
since each of two subsets is partitioned with respect to the same dimension 
or feature, whereas in general two different subsets would be partitioned with 
respect to different dimensions. The variables can easily be plotted with 

respect to an orthogonal coordinate system that displays the major clusters. 

O 
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Ill, Relations Among Linear Contrasts Based on Different Subsets of the Stimuli 
Suppose that a set of n stimuli hats been ^partitioned into two subsets, 
and Sg , according to a "contrast component" b derived by Method A or 
Method B. Now , which contains, say, n^ stimuli, can be further sub- 
divided according to a contrast component b^ . 

An important question, first raised by McQuitty (1967 ^ cf » > his Figure l), 
is, "How are b and b^ related?" In particular, are b and b^ orthogonal 
or not? At first this question seems meaningless, since b has n elements, 
while b^ has n^ < n elements. Ho\^ can a scalar product be computed between 
two vectors that have different numbers of elements? 

The following solution to this problem takes advantage of the fact that 
the stimuli have been partitioned with respect to identifiable underlying linear 
contrasts . 



Let R = 




Ri 2 




r— 1 

CVJ 


R 22 



be the partitioned similarities matrix. The vector 



b has been computed using R , while the vector b^ is based only on R -) 1 . 
A vector d , of n elements, can be constructed such that the first n^ 
elements of d are proportional to the elements of b^ , and d'Rd is maxi- 
mized under the further constraint d’d = 1 . 

Let d' = , where d^ = kb^ for some unknown scalar k . . The 

problem is to find dg and k such that 



[kb| d^] 

(27) F(k,dg) = 

is maximized \inder the constraint k“" + 
method of Lagrange inultipliers , let 




Ri2 






CVJ* 

CVJ 

p:: 






^^2 = 


= 1 . 



Once again applying the 
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( 28 ) G(k,d2,7v) = ck^ + d^R22l2 -^2 " ' 

where c = b^'R^ .| b.| and v = ^21-1 

Differentiating G with respect to k , , and X , and setting the 

derivatives equal to 0 yields 

(29) = 2R2gd2 + 21^ - 2Xdg = 0 ; 

(50) gg = 2ck + Sv'dg - 2Xk = 0 ; 

(51) g^ = 1 - k^ - d^2 = 0 . 

It may be safely assumed that k / 0 . Premultiplying (29) by ^ and 
dividing by 2 yields 

(52) ^2^22^2 ^2- ■ ^"^2 " ° * 

Multiplying (50) by k and dividing by 2 yields 

(55) ck^ + lor'dg - = 0 . 

Adding (52) and (55) yields 

(5U) ck^ + 1^R22^2 ^ -2 " -^2^ " ° • 



(55) 



(56) 



Equations (5^) and (5l) imply 

X = ck^ + ^^^22^ 2kv'dg = F(k,d2) 

Now ( 29 ) and (50 ) cein be rewritten more conveniently as 



I >22 Z 




'^2 




r n 

i2 








= X 




V* c 




k 




k 


- ^ 
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From (51), (55), and (56), it follows that 



is the normalized latent 



vector corresponding to the largest latent root 

l>nien the values of dg and k have been 
of accuracy, then the solution to the original 

(57) 1' = ^2^ > where d^ = kb^ 

Now d , unlike b^ , has the same number 
of the angle betv/een b and d is simply ^*R^[(b’Rb)(d'Rd)]2 . 

IV. Example 

The similarities matrix for this example (Table l) contains phi - coefficients 
between pairs of subjects, based on 90 binary responses. The data were obtained 
in a study of reading comprehension. There were three "treatments,” indicated in 
the tables and figures by X , Y , and Z . One of the basic hypotheses of the 
study was that subjects within a treatment group vfould be relatively homogeneous 
in terms of their response patterns, and that the groups would be distinct. In 
other words, there was an a priori three-cluster hypothesis with respect to the 
subject -space. 



^2 



of 



R22 i 



determined to the desired degree 
problem is 



of elements as b . The cosine 



Insert Table 1 about here 

Table 2 gives the coordinates of the subjects on four contrast components, 
namely the first and second Method A and the first and second Method B com- 
ponents (Al, A2, Bl, and B2, respectively). In each case, the subjects have 
been ordered with respect to their coordinates. Differences between suc- 
cessive pairs of coordinates are given, and partitions have been made according 
to the "weak criterion" suggested in Section II. 



The right-hand column of Table 2 indicates the first partition of the sub- 
jects according to IICA. Note that this split is the only one which fails to 
conform to the initial hypothesis. Furthermore, this mis classification by IICA 
is rather "robust," persisting even when the similarities are converted from 
phi- coefficients to rank scores . 

Insert Table 2 about here 

The extension of A2 and B2, from a 19-subject subspace to the full 28- 
subject space, is outlined in Table 5. k was found to be 0.2675 for Method A 
and 0.2705 for Method B. The cosine of the angle between the first and extended 
second components was found to be -.7755 for Method A and O. 786 O for Method B. 

Insert Table 5 about here 

Figure 1 shows a plot of the subjects in the plane determined by the first 
and extended second components according to Method A. Figure 2 shows a similar 
plot for the Method B components . The appropriate columns of Table 5 have been 
scaled according to b'Rb or d'Rd . 

Insert Figures 1 and 2 about here 

Figures 5 and 4 show plots of two-dimensional orthogonal solutions obtained 
by Methods A and B, respnctively. The three-cluster structure is apparent. 

Insert Figures 5 and 4 about here 
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V. Summary and Conclusions 

McQviitty's (1967) technique of hierarchical cluster analysis --Iterative 
Intercolumnar Correlational Analysis — has been discvissed in terms of matrix 
algebra and geometry. Under this interpretation, it has been shown that IICA 
achieves a solution by transforming the stinralus -vectors themselves toward a 
bipolar, one-factor struct\ire. Two alternative methods were suggested for 
extracting a single bipolar factor directlj.' from the initial similarities 
matrix. Extension of linear contrasts from smaller to larger subspaces was 
also discussed. The major features of the new methods were illustrated in the 
analysis of some data from a reading comprehension study. 
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Table Headings 

Table 1. Similarities matrix. Hii -coefficients between pairs of subjects, 
based on 90 binary responses. The partitioning indicates the a priori three- 
cluster hypothesis. (Decimal points are suppressed. ) 

Table 2. Coordinates of subjects on contrast components. Subjects have 
been ordered with respect to their coordinates, and differences between succes- 
sive pairs of coordinates are given in the colximn to the right of the coordinates 
themselves. Partitions according to the weak criterion are given. The right-hand 
column gives the first partition according to IICA. (Decimal points are sup- 
pressed. ) 

Table 5. Extension of second Method A and Method B contrast components. 
Blanks in columns A2 and B2 are filled by the extension ^ . Other elements in 
columns A2E and B2E are equal to k times the corresponding elements of columns 
A2 and B2, respectively. 
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O rH-d--d-C0 
O ITN ITN ITN -d- ITNVO 


iTN-d- ITN ITN Q rH 00 C0» Q 
CVJCVJCVJrHCO 1 iHrHKNCVJ 


ON 

X 


CVJ C^CO CVJ ON O O 

CO iTN-d* -d- VO lOv Q 


^ Q t^co g> t^oo H w 
to to CVI CVI H H 


O ITN ITN-d- KN -d- ON KNVO 

CVJ H CVJ CVJ CVJ H 




H KNK^VO C^CVJ Q O 
^ VO VO VO VO ITVVO O ^ 


lOO-d-lOCVIOHO\Og\CM ^-^ovcy ovio 

CVJrHrHrHKN CVJrHKNrHCVJrHrHrH 1 CVJrHrH 


u 


«d*VOCVJCVJVOi/NOCVJ 
^ VO VO ITN C^vo o VO -d* 


^ CVI H VO CO (TWO CO to VO to cy ^ to Q to VO lO CVI 


VO 

X 


VOCOCO-d* C^OI^C^ON 
ITNVO VO KNCO O VO LTN LTN 
rH 


CVI lotovo -=f Hioovcvi ONtoco loco Q CTvco cy ^ 

H H H 1 H H 1 r-t 

1 ' 


if\ 

X 


ON VO VO O C^vo VO CVJ 

ITN c^-d* C CO C^vo VO 


K^toHt^o^l-c^cyHcyc^o-d■cywoMocyt^ 

HHcycy iHcyH 




CT\ H O C^-d* CVJ KNCO 
-d* VO O -d* KNITNVO -=f 
rH 


cy On to lO J- o to to«to C3\ on to on to vo H On to 
HcycyHtoHtoHcy h h ih • 


X 


COVOO HVOCOCVJKNt^- 
ITN O VO C^VO VO VO -=t- 

rH 


-cjnc^ioco loocyHpo cyiocyH^ c^-d- h vo 
H H cy cy H cytocy h h 




-=f ovo C^VOCOVO H t> 
ITN o C^-d* C^VOVO VO IT 
r=5 


- CO tovo f^ONC^t-OJpCO 

xHHHcyH icycycy HHiHcynHcynn 


H 

X 


o -d-CO ON ONVO-d* 0^ 

O KN i/N-d* KN i/N-d* -d* CC 

r=H 


1 ON ONVO tovoHtopcyHQQ oJOJiogc^cy^ 

I toj- toitocycyitocycycy h h 




o 

iH CVJ KN ji± liNVO h-“CO ON H CVJ KN-d* liNVO ^OD PN H CU iO\Q 



w 



V 
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Table 1 (Continued) 





Z1 


Z2 


Z5 


Z4 


Z5 


Z6 


Z7 


Z8 


Z9 


ZIO 


XI 


-1 


50 


20 


22 


22 


5 


10 


7 


2 


15 


X2 


8 


17 


l4 


U 


22 


-12 


15 


21 


12 


17 


X5 


0 


12 


5 


2 


1 


-15 


7 


4 


1 


6 


X4 


5 


19 


9 


15 


9 


-5 


l6 


1 


9 


-5 


X5 


2 


7 


0 


4 


2 


-12 


9 


5 


2 


7 


X6 


9 


15 


8 


5 


8 


-10 


9 


8 


-2 


15 


X7 


-6 


-5 


-2 


-4 


-5 


-20 


5 


-6 


-5 


2 


X8 


19 


22 


17 


14 


19 


-2 


22 


l4 


19 • 


5 


X9 


0 


25 


15 


24 


25 


7 


24 


9 


5 


16 


Yl 


25 


24 


25 


15 


20 


-1 


18 


18 


50 


24 


Y2 


7 


29 


25 


11 


19 


7 


5 


11 


19 


21 


Y5 


-7 


18 


15 


10 


19 


6 


-2 


-5 


-1 


12 


y4 


-1 


5 


-5 


-19 


-20 


-l4 


-12 


-5 


0 


-4 


Y5 


20 


57 


48 


50 


56 


27 


28 


26 


56 


40 


y6 


-5 


12 


9 


-8 


-7 


5 


-1 


-16 


-7 


-1 


Y7 


15 


40 


29 


18 


18 


15 


15 


25 


28 


22 


y8 


17 


11 


18 


. 14 


15 


-9 


10 


19 


25 


18 


Y9 


20 


55 


28 


18 


17 


8 


20 


15 


2T ■ 


21 


Z1 


100 


58 


52 


49 


45 


19 


58 


71 


59 


40 


Z2 


58 


100 


55 


29 


26 


25 


56 


56 


26 


29 


Z5 


52 


55 


100 


45 


65 


25 


47 


55 


50 


58 


Z4 


49 


29 


45 


100 


48 


12 


61 


51 


55 


42 


Z5 


45 


26 


65 


48 


100 


51 


4l 


57 


58 


71 


Z6 


19 


25 


25 


12 


51 


100 


24 


18 


4 


27 


Z7 


58 


56 


47 


6i 


4l 


24 


100 


55 


4i 


48 


Z8 


71 


56 


55 


51 


57 


l8 


55 


100 


75 


50 


Z9 


59 


26 


50 


55 


58 


4 


4i 


75 


100 


52 


ZIO 


40 


29 


58 


42 


' 71 


27 


48 


50 


52 


100 



ERIC 
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Table 2 

Coordinates of Subjects on Contrast Components 



O 

ERIC 



S# 


P2 . 


diff 


S# 


A 2 


diff 


Z 1 


-259 




y 4 


-525 








005 






055 


Z8 


-254 




Y 5 


-288 








000 






015 


Z 9 


-254 




y6 


-274 








010 






018 


Z 5 


-225 




Y 2 


-256 








001 






006 


Z 5 


-222 




Y8 


-250 








008 






051 


ZIO 


- 2 l 4 




Y 9 


-219 








026 






001 


Z 7 


-188 




Y 1 


-218 








012 






027 


Z 4 


-176 




77 


-191 








005 






126 


Z 6 


-171 




Y 5 


-066 








055 






154 


Z 2 


-116 




Z 2 


069 








010 






056 


Y 5 


-106 




Z 6 


124 








047 






077 


Y 1 


-059 




Z 5 


202 








001 






Oil 


Y 9 


-058 




Z 9 


215 








012 






005 


Y 7 


-047 




ZIO 


215 








020 






002 


Y 2 


-027 




Z 4 


217 








017 






057 


Y 6 


-010 




Z 5 


254 








009 






005 


y8 


-002 




Z 7 


257 








048 






002 


Y 5 


046 




Z 1 


259 








064 


■ 




017 


y 4 


110 




Z8 


276 








080 








X 9 


190 














000 








XI 


190 














001 








X 4 


191 














005 








X8 


19 *<- 














04 l 








X 2 


255 














054 








x 6 


270 














016 








X 5 


286 














018 








X 7 


505 














006 








X 5 


509 











S# 



B 1 



diff 



X 5 -512 
X 7 -505 
X 5 -288 
x6 -271 
X 2 -258 

X 8 -197 

Xl 4 - -195 
XI -195 
X 9 -195 
-112 
Y 5 -050 

Y 8 -002 



007 

016 

017 
055 
04 i 
oou 
000 
001 
081 
062 
048 







012 


y 6 


009 


015 


Y 2 


025 


019 


Y 7 


o 42 


012 


Y 9 


054 


002 


Y 1 


055 


047 


Y 5 


102 


012 


Z 2 


115 


058 


Z 6 


172 


005 


Z 4 


175 


012 


Z 7 


186 


025 


ZIO 


212 


008 


Z 5 


220 


001 


Z 5 


221 


012 


Z 9 


251 


001 


Z 8 


252 


006 


Zl 


258 





s# 


B 2 


diff 


IICA 


y 4 


-525 


054 


Z 6 


Y 5 


-289 


014 


Z 2 


y6 


-274 


017 


Z 9 


Y 2 


-257 


006 


Z8 


Y8 


-251 


051 


ZIO 


Y 9 


-221 


001 


Z 5 


Y 1 


-220 


027 


Z 5 


Y 7 


-195 


125 


Z 7 


_Y 5 __ 


_-o67 




Z 4 






155 




Z 2 


068 


057 


Zl 


z6 


124 


C76 


Y 7 


Z 3 


200 


Oil 




Z 9 


211 


005 


Y 9 


ZIO 


214 


005 


_Y 1 _ 




217 


056 


'Ve 


Z 5 


255 


005 


y 4 


Z 7 


256 


002 


y8 


Zl 


258 


017 


Y 5 


z 8 


275 




Y 2 



X 8 

X4 

X 7 

X 5 

x 6 

X 5 

X 2 

XI 

X 9 
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Table 5 

Extension of Second Method A and Method B Contrast Components 



S# 


A 1 


A 2 


' A 2 E 


B 1 


B 2 


B 2 E 


XI 


190 


— 


-302 


-193 


— 


-302 


X 2 


255 


— 


-332 


-238 


— 


-332 


X 5 


286 


— 


-349 


-288 


— 


-349 


X 4 


191 


— 


-279 


-193 


— 


-279 


X 5 


310 


— 


-365 


. -312 


— 


-564 


X 6 


270 


— 


-322 


-271 


— 


-322 


X 7 


303 


— 


-329 


-305 


— 


-528 


X 8 


194 


— 


-307 


-197 


— 


-306 


X 9 


190 


— 


-298 


-193 


— 


-298 


Y 1 


-060 


-219 


-059 


055 


-220 


-060 


Y 2 


-027 


-256 


-068 


025 


-257 


-070 


Y 5 


046 


-288 


-077 


-050 


-289 


-078 


y 4 


no 


-325 


-086 


-112 


-323 


-087 


Y 5 


-106 


-066 


-018 


102 


-067 


-018 


y 6 


-010 


-274 


-073 


009 


-274 


-074 


Y 7 


-047 


-191 


-051 


042 


-193 


-052 


y 3 


-002 


-250 


-067 


-002 


-251 


-068 


Y 9 


-058 


-219 


-059 


054 


-221 


-060 


Z 1 


-239 


259 


069 


258 


258 


070 


Z 2 


-116 


069 


018 


113 


068 


018 


Z 5 


-223 


201 


054 


221 


200 


054 


Z 4 


-176 


217 


058 


175 


217 


059 


Z 5 


-222 


254 


068 


220 


253 


068 


Z 6 


-171 


124 


033 


172 


124 


034 


Z 7 


-188 


257 


069 


186 


256 


069 


Z 8 


-234 


276 


074 


232 


275 


074 


Z 9 


-234 


213 


057 


231 


211 


057 


ZIO 


-214 


215 


058 


212 


214 


057 
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Figure Captions 

Fig. 1. Plot of columns A1 (b) and A2E (d) from Table 5. The 

coordinates have been multiplied by the square root of ’'variance -account ed-f or. ' 

The cosine of the angle between the axes is -.76» 

Fig. 2. Plot of columns B1 (b) and B2E (d) from Table 3. The coordi- 

nates have been multiplied by the square root of "variance -accounted-f or. " The 

cosine of the angle between the axes is 0.79* 

Fig. 5. Two-dimensional orthogonal solution, Method A. 

Fig. 4. Two-dimensional orthogonal solution, Method B. 
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